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Long-range frustrating interactions, even if their strength is infinitesimal, can give rise to a dramatic 
proliferation of ground or near-ground states. As a consequence, the ordering temperature can 
exhibit a discontinuous drop as a function of the frustration. A simple model of a doped Mott 
insulator, where the short-range tendency of holes to phase separate competes with long-range 
Coulomb effects, exhibits this "avoided critical" behavior. This model may serve as a paradigm for 
many other systems. 



A wide variety of systems display equilibrium domain 
patterns characterized by periodic (or nearly periodic) 
variations of an order parameter. These patterns are 
stabilized by competing interactions. Linear arrays of 
stripes and hexagonal arrays of bubbles are ubiquitous 
in thin films of magnetic garnets, and ferrofluids. Simi- 
lar morphologies are also seen in Langmuir films, mem- 
branes, semiconductor surfaces, and many other systems 
in which an otherwise uniform ground state is thwarted 
by a competing "frustrating" interaction of one sort or 
another^. Lately, stripe structures have been detected 
I 2 1 in doped Mott insulators, including the high T c su- 
perconductors: the ordered states in these compounds 
consist of arrays of charged stripes which form antiphase 
domain walls between antiferromagnetically ordered spin 
domains. In the absence of a frustrating Coulomb in- 
teraction (i.e. for neutral holes), a lightly doped Mott 
insulator is unstable to phase separation into a hole- 
rich "metallic" phase and a hole-deficient antiferromag- 
netic phase. Electrostatic repulsions forbid macroscopic 
charge seperation; the compromise leads to the formation 
of stripe morphologies on an intermediate scale. 

In this letter, we shall explore the effect of fluctua- 
tions on the periodic structures in a simple model of uni- 
formly frustrated 0(n) spins. We study the problem us- 
ing two complementary approaches: a low temperature, 
spin-wave expansion, and a perturbative expansion for 
large n, which we carry through to order 1/n 2 . It is 
found that the ordering temperature, T C (Q), as a func- 
tion of the strength of the frustrating interaction "Q" 
may, in certain instances, satisfy the inequality 



T c {Q = 0) >Hm T C (Q). 



(1) 



In other words, an infinitesimal amount of frustra- 
tion depresses the ordering temperature discontinuously! 
Specifically, we shall argue on the basis of a low temper- 
ature expansion about stripe like ground states, that, in 
the absence of lattice effects, T C (Q) — for n > 2 and 
Q > 0. Lattice anisotropies elevate T C (Q) from zero, but 
the discontinuity in T C (Q) persists for 2 < d < 3 and 
n > 2. For n = 2, the lower critical dimension is three, 
and here the finite temperature "ordered" phase exhibits 
power law decay of correlations; the model with n = 2 



has the same hydrodynamic description as a smectic liq- 
uid crystal. 

For low frustration Q, the behavior of the system is 
controlled by the proximity to the "avoided critical tem- 
perature" T C (Q — 0). The following picture emerges of 
the thermal evolution of the model, as summarized in 
Fig. 1: At temperatures somewhat above T C (Q = 0), 
two large lengths govern the exponential decay of corre- 
lations. As the temperature is lowered below a crossover 
temperature T\(Q) ~ T C (Q = 0), the system enters a low 
temperature regime characterized by an oscillatory spin 
structure function with a single length controlling the ex- 
ponential decay of correlations at long distances. T\(Q) 
is a "disorder line" in the sense that as T approaches Ti 
from below, the wave-length of the oscillations diverges, 
but no phase transition occurs. As T is lowered fur- 
ther, the wave-length decreases until, as T — > T C (Q), it 
smoothly approaches the period of the ordered phase that 
appears below T C (Q). However, an additional crossover 
occurs at a temperature T2(Q) which lies between T C (Q) 
and Ti(Q), such that for T 2 (Q) > T > T C (Q), there are 
again two long lengths characterizing the fall-off of cor- 
relations, where the new length is akin to the Josepshon 
length in the ordered phase of the unfrustrated system. 
This second length can only be seen in the context of a 
1 /n expansion. The existence of multiple correlation and 
modulation lengths is a common feature of the physics 
of all the various frustrated systems alluded to above M . 
A finite size scaling analysis, which is a simple extention 
of an argument presented previously' 9 ! allows us to iden- 
tify the longest length in the temperature regime above 
T\ and below T 2 as a "domain" size, R, within which 
the physics is essentially that of the unfrustrated system, 
and to extract the scaling relation (which is reproduced 
by the large n results) 



(2) 



where £o ~ [T C (Q = 0) — T\~ v is the correlation length 
in the unfrustrated system at temperature T. 
The Coulomb Frustrated Ferromagnet: As a con- 
crete example, we consider a system with a short-range 
tendency to phase separation which is frustrated by a 



1 



long-range Coulomb interaction. A simple spin Hamilto- 
nian which represents these competing interactions is 

Q ^ S(x)S(y) 
2 



Hn = - 



E 

<x,y> 



S(x)S(y) 



E 



y\ 



(3) 



Here, S(x) is a coarse grained scalar variable which rep- 
resents the local charge density. Each site x lies on a 
cubic lattice (of size N) and represents a small region 
of space in which S(x) > 0, and S(x) < correspond to 
the positively and netatively charged phases respectively. 
The first "ferromagnetic" term represents the short-range 
(nearest-neighbour) tendency to phase-separation, while 
the second term is the Coulomb interaction. Non-linear 
terms in the full Hamiltonian typically fix the locally 
preferred values of S(x). One may consider d ^ 3 di- 
mensional variants wherein the spins lie on a hypercu- 
bic lattice, and the Coulomb kernel in Hq is replaced by 
Q\x — y\ 2 ~ d ■ H can be fourier transformed as 



H = j2m\s(k)f 



where the kernal 



J(fc) = -[Qk- 2 + r + k 2 



(4) 



(5) 



where tq = —2d, and the ellipsis represents higher order 
terms in powers of k. We will neglect these terms for 
now, as they are unimportant in the contiuum; however, 
we will need to include some of these terms when we 
treat lattice effects since they are the ones that reduce 
the full rotational symmetry of free space to the point 
group symmetry of the lattice. 0] 
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FIG. 1. Schematic Phase Diagram. Q is the strength of 
the frustration. For T > 7i(Q) and T 2 {Q) > T > T C (Q), 
there are two long lengths governing the fall of correlations. 
For T < Ti(Q), the spin-spin correlation funcition exhibits 
long distance oscillatory structure. The thick black dot marks 
T C (Q = 0), the ordering temperature in the absence of frus- 
tration; this is what we term "the avoided critical point". In 
the continuum limit, T C (Q > 0) = 0; lattice effects result in a 
non-zero, yet small T C (Q). 



We now generalize this model, allowing the spins to 
have n components, and replacing all two spin products 
in Hq with a scalar product. We treat both the "soft- 
spin" version of this model, in which we include the non- 
linear interaction 



H 



soft 



; £[S 2 (oQ - If 



(6) 



with u > 0, or the "hard-spin" version, which can be 
viewed as the u — > oo limit of the soft-spin model, in 
which we instead enforce the local constraint, |S(a?) | = 1. 

When n > 2, we can construct a set of ground-state 
configurations which are simple spirals of the form 



S 9 (x) — acos(fc min • x) + bsin(fc min • x), 



(7) 



where k m i n denotes any which wave-vector which mini- 
mizes the interaction kernel, J(k), in Hq, and the pref- 
actors satisfy 



a-b = 0; aa=bb = l. 



(8) 



It is readily seen that such states are unstable to trans- 
verse fluctations. One may expand H so f t to quadratic 
order in fluctuations about the ground state AS(ir) = 
[S(x) — S 9 (x)], and estimate the thermal average of 
AS 2 (if). For d — 3 and a vanishing lower cutoff e on 

I k | k m in 1 1 : 



< (AS) 2 > ( n - 2)VQ 1 
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Q 1/4 ln| 



(9) 



For n > 2, the leading order divergence is 0(e _1 ), in- 
depepent of d; for XY spins (n = 2) the leading order 
divergence, in 2 < d < 3, is 0(e d ~ 3 ). In d = 3, XY spins 
exhibit power law correlations at low temperatures. As 
we shall show, for small Q and n — 2, 3, these simple spi- 
rals are the only possible ground states; for n > 4 other 
types of "multi-spiral" ground-states are possible. 
The Three-Dimensional Spherical Model: To make 
the phase diagram non-trivial, yet tractable, we may 
solve the scalar spin model subject to the single mean 
spherical constraint! 1 ' 5 ] 



E < * 2 ( f ) >= 



N. 



(10) 



It is known [11] that, in many respects, the spherical 
model is equivalent to the n — > oo limit of the 0(n) 
model. Here, the effective Hamiltonian is the same Hq 
defined above, with tq — > r, where r is a Lagrange multi- 
plier determined implicitly from the constraint equation 
( |l0| ) . In order that all modes have a bounded Boltzmann 
weight it is necessary that r > —2\[Q. By equipartition, 



< \S(k)\ 2 >= T[k 2 + Q/k 2 + r] 1 , so the mean spherical 
constraint reads 



1 

T 



\k\<A 



d 3 k 

(2tt) 3 fc 5 " 



Qk- 



(11) 
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where A is an ultraviolet cutoff. If this equation cannot 
be satisfied for any value of r > — 2 V / Q, then the system 
is at or below criticality. 

Observe, from Eq.(|ll|), that for Q > 0, the integral 
diverges when r — ► —2\fQ. Thus, the constraint can be 
satisfied for any non-zero T: T C (Q > 0) = 0. By contrast, 
when Q = ( which is the standard three-dimensional 
short-range ferromagnet) T c is non-zero. Thus a discon- 
tinuity in T C (Q) is seen to exist. However, even though 
T c = for Q > 0, there is a genuine zero temperature 
phase transition with the usual n — > oo critical expo- 
nents, e.g. v = 1 and 7 = 2 in d = 3. As we shall see, 
lattice effects elevate T C (Q), but do not change the criti- 
cal properties, nor, in d = 3, eliminate the discontinuity. 

The pair correlator is given by 



G(x) 



1 



(2tt) 

T 



., , d 3 k < \S(k)\ 2 > exp[ik-x] 



dk,^ In *^L (12) 



27r 2 \x\J (k 2 + a 2 )(k 2 + /3 2 ) 



where 



(13) 



When r > 2\fQ, the integral can be readily evaluated by 
applying the residue theorem to the poles lying on the 
imaginary axis at k = ±ia, ±i/3, 



G(x) 



T((3 2 e 



a 2 e -a|x| 



4ir\x\((3 2 ~ a 2 ) 



(14) 



Note the existence of two macroscopic correlation lengths 
- a consequence of charge neutrality: In Hq, the spins 
portray charges, and therefore must sum to zero, 



G(x 



d 3 3 



-< 



S(x] 



d 3 x\ 



>=0. 



(15) 



Whenever G is dominated by its long-distance behavior, 
the integral can vanish only if G(x) contains positive and 
negative contributions, as in Eq.dl4). The latter integral 
can be made to vanish only if G(irjcontains, at least, two 
length scales. At high temperatures, the length 



\Re{(3}\- 



-- 1/2 { for r > 2y/Q). 



(16) 



plays the role of the correlation length of the canonical 
short-range ferromagnet (i.e. with Q — 0). Note that 
now, however, an additional correlation length appears: 



6 = \Reia}]- 1 « Q~ 1/2 /£ 



(17) 



Thus, ^2 i n the limit of weak frustration, Q <C 1. 

The analytic continuation of Eq.(|lJ) to low tempera- 
tures, r < 2y/Q, is 



T 

G(x) = — exp(— ai|x|) 

47T 



(o 



i sina2|a?| + 2aitt2 cosc^!^ 



4aia2 x 



(18) 



where a = ol\ + z«2- The temperature T\, defined by 
r(T = T\) = 2^/Q, marks a dramatic crossover. At 
low temperatures (T < Xi), the system possesses a sin- 

gle correlation length £ = = 2[r + 2^/Q] , 

and a single modulation length Ld = 27r/|ct2 1 = 



4tt[ 



-1/2. 



at high temperatures (T > Ti), the 



system possesses two distinct correlation lengths. When 
T = , the modulation length diverges as ~ 
(Ti — T)~ x / 2 . Many quantities of interest (e.g. the 
specific heat Cy) albeit analytic, display a crossover at 
T\(Q). As Q tends to zero, the crossover temperature 
T\(Q) tends to T C (Q = 0). Thus despite the nonex- 
istance, for Q = + , of a phase transition at or near 
T C (Q = 0), the system is governed, in part, by the prox- 
imity to the avoided critical temperature T C (Q = 0). 
Avoided Critical Behaviour To 0(n 2 ) : We will 
now examine corrections to the spherical limit. In a (1/n.) 
expansion^ 4 ! , the soft term of constraint, [H so f t — Ho) is 
taken to be small with u = 0(1 /n) > 0. The pertur- 
bation theory in u is then selectively resummed treating 
1/n as the small parameter. As our unperturbed Hamil- 
tonian, we take Hq in Eq. (3) with a temperature depen- 
dent chemical potential, ro + 2^/Q, which changes sign 
at T — Tmf, and is increasingly negative at low T. The 
Dyson equation implies 



G" 1 (fc) = Go 1 (fc) + S(fc) 



(19) 



where, in the continuum limit, G 1 = ro+k 2 +Qk 2 , and 

[—£(&)] is the self-energy. T c , if it exists, is determined 
implicitly from the solution of the equation 



mmiG-^k)} = 0. 

k 

To zeroth order in 1/n: 

= r + k 2 + Qk- 2 . 



(20) 



(21) 



Here r — tq + E°, where [— S°] denotes the self- 
consistently computed O(n ) correction to the self- 
energy. At low temperatures, 



E°(Q, r) « E°(Q = 0, r = 0) + — 



7T V r + 2VQ' 



(22) 



where E°(Q = 0,r = 0) is the value of the O(n ) self- 
energy at criticality for the standard three-dimensional 
ferromagnet and the second term is the 0(n ) self-energy 
of a one-dimensional spin chain with a nearest neighbour 
exchange interaction proportional to 1/Q. Since E° is 
manifestly positive, and diverges as r — > r m ; n = —2yfQ, 
Eq. ( pp[ ) is satisfied only in the limit r = —00; to this 
order T C (Q > 0) = 0. We have extended this analy- 
sis H to 0(n~ 2 ). By evaluting diagrams self-consistenly, 
one observes that all 0(n _1 ) and 0(n~ 2 ) self-energy con- 
tributions are explicitly positive, or cancel against more 
divergent positive contributions. We outline here how 



3 



this is done to 0(n x ). In Fig. 2, the fc-independent 
E° > is the single zeroth order (O(n )) contribution. 
To 0(n _1 ) there are two additional diagrams: T, A (k) 
and the k- independent E s . Inserting, self-consistently, 
G(k) = [Gv\k) + E° + E A (fc)]- x into the integral ex- 
pression E° = J pjj^ G(k), automatically generates E B , 
as well as higher order diagrams. This integral diverges 
as r — > r m i n . The self-energy E is positive, and thus 
can only further thwart any tendency to order. A similar 
anaylsis [5] may be repeated to 0(n~ 2 ). 
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FIG. 2. The Self Energy Corrections. The thin dashed 
lines denote bare interaction, the thick dashed lines represent 
dressed interactions (i.e. bare interactions screened by a geo- 
metric series of bubble diagrams) , and the solid lines denote 
propagators. 

All this indicates that to 0(n~ 2 ), r = r m ; n is attain- 
able only at T c = 0. Unfortunately, the results of the 1/n 
expansion are not independent of Q for small Q; we can- 
not safely draw conclussions concerning the Q — ► limit 
as there could, in principle, be a change in the behavior 
of the system when Q ~ 1/n. Nevertheless, the results 
strongly support the contention that the n — > 00 limit is 
not singular, and that the spherical model captures the 
important physics of the system for any large n. 
Algebraic Crossover: We have also computed the pair 
correlator to (^(n" 1 ). At very low temperatures, where 
< r + 2y/Q *C y/Q, the propagator at intermediate 
momenta (1 » |fc| > Q 1/A ) is 

G-^fc) w k 2 + [2Y?/(n 2 u)]\k\ + r + Qk' 2 . (23) 
We note that 



G(\x\) ~ |fr 2 for £j < |f| < L D < £ 
G(\x\) ~ If!" 1 for |f| « £ J; 



(24) 



where the correlation length £j = n 2 u/(2E°), is defined 
in a way analagous to the Josephson length l 6 l in the or- 
dered phase of a system with Goldstone modes. Thus 
at sufficiently low temperatures, T < T2, the (nonoscilla- 
tory) spatial behavior of the correlators is again governed 
by two length scales. 

Lattice Effects: The fact that the lattice system has 
discrete, rather than continuous rotational symmetry is 
reflected in higher order terms in powers of k, in the ker- 
nal J(k) in Eq. (||); the lowest order term of this sort 

previously |3J for n 



but the avoided critical phenomena, i.e. the fact that 
limg^o T C (Q) < T c (0), survives for 2 < d < 3. More 
generally, if we apply the above spin-wave analysis and 
a Lindemann criterion for T c , then the same calculation 
leads to the conclusion that, once again, the large n re- 
sults are qualitatively correct for finite n. 
Multi-spiral states Whenever n > 2, any ground state 
configuration can be decomposed into Fourier compo- 
nents, S 9 (x) = J2tLi { a * COf #mL ■x]+h, sin[fc^ n ■ x]j 

where k l min are chosen from the set of wave vectors which 
minimize J(k). So long as these wave- vectors are "non- 
degenerate" , in the sense that the sum of any pair of 
wave vectors, k l min ± k J min is not equal to the sum of any 
other pair of wave vectors, and "incommensurate" in the 
sense that for all i and j, 2(k l min + k 3 min ) is not equal to 
a reciprocal lattice vector, it is straightforward to prove 
H that the condition [S s (af)] 2 = 1 can be satisfied only 
if M < n/2. (These conditions are always satisfied for 
Q < 4.) Thus, for n < 3 only simple spiral [M = I) 
ground-states are permitted, while for n — 4, a double 
spiral saturates the bound ^ 12 1 . 
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00; they produce a T C (Q > 0) > 0, 
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